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ABSTRACT 

By use of distilled water samples in the laboratory, and very clear lakes 
in the field, a technique has been developed where the atmosphere and sur- 
face noise effects on LANDSAT signals of water bodies can be removed. The 
residual signal is dependent only on the material in the water. . The gen- 
eral type and concentration of the material can be determined from LANDSAT 
tapes. When this material in the water is living algae or weeds, its con- 
centration is related to the enrichment or eutrophication of the lake. 

The algae and weed biomass is at a maximum in late August which is there- 
fore the optimum time to categ'brize lakes. 

For this project, the Wisconsin Department of Natural Resources helped pro- 
vide water samples, and the University of Wisconsin under NASA support helped 
provide laboratory and computer analysis and aerial observation of the 
lakes. The Bendix Multispectral Data Analysis System (M-DAS) provided 
a Color Categorized Image, of several hundred lakes. These lakes were 
categorized for tannin or non-tannin waters and for the degrees of algae, 
silt, weeds, and bottom effects present. The initial categorization is 
being field checked in late August 1975. 
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INTRODUCTION 


When LANDSAT images are used for mapping land targets essentially two fac- 
tors are involved: (1) How the sun and skylight interact with the target 
to form a reflected signal, and (2) how the atmosphere effects this signal 
before it reaches the sensor. Land targets are considered as diffuse re- 
flectors and Lambert's Law for diffused reflection can be used in the anal- 
ysis. 

IVhen LANDSAT images are used for mapping or monitoring water quality the 
process becomes m.ore complex. Sunlight and skylight interact with water 

surface, the particles in the water volume , _and in som encase s_with_ the 

bottom . The signal from the water surface must be treated as primarily 
specular reflection but there is a small surface diffused reflection com- 
ponent caused by dust, foam and other contaminants on the water surface. 

The surface specular component varies greatly depending on the skylight 
condition and is essentially independent of the material in the water 
volume (1). The signal from the water volume can be treated as diffused 
reflection from the suspended particles in the water. The signal from the 
bottcm is troublesome noise which follows the laws of diffused reflection. 

All of these signals combine to form the total signal from the lake or other 
water body. This total signal is modified by the atmosphere prior to its 
reaching the satellite sensor. 

The most important component of the signal from the water is that caused 
by the material in the water volume. This is called volume reflectance 
or backscatter and is herein denoted by p,. . There is a different p^- for 
each color (wavelength) . The volume reflectance for a water is primarily 
caused by light being diffusely reflected from suspended particles in the 
water between the water surface and the depth where the energy is extinguish- 
ed. Pure or distilled water is assumed to contain no suspended material? 
therefore the volume reflectance of the suspended particles in distilled 
water is essentially zero. Ifowever, even in the laboratory there will be 
some diffuse backscatter from distilled water caused by the water molecules, 
dust, foam, or other impurities on the ivater surface. Let the laboratory 
diffused reflectance from impurities on the water surface be indicated by 
Pj.^. There is a different pg^^ for each color or wavelength of energy. 

As material is added to a pure ivater sample the only factor that will be 
significantly altered will be the volume reflectance, . The dust, foam 
and other impurities arelikely to remain unchanged so the pcj should not 
change. In some cases oil slicks can alter this psl but in all other cases 
it is considered unchanging. It is p^, volume reflectance, which changes 
as material is added to the water. The type and amount of material added 
to pure water is what alters the water quality. Each type of material such 
as red clay, green taconite mine tailings, blue-green algae, etc. reflect 
differently at different wavelengths. These different materials have unique 
spectral reflectance signatures which are indicated by the Py at different 
wavelengths. The LAJ®SAT Satellite has four different sensor bands so the 
type of material in water should be detectable frcm the four bands of LANDSAT 

For a particular size and shape of particle, as more material is added to 
the water there vdll be more particles to backscatter the light. The vol- 
ume reflectance, P^, will increase as will the total signal sensed by the 
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Satellite. Therefore, for a particular type of material, its concentration 
should be related to the LANDSAT signal strength. 

It should be ]x)ssible to determine both the type of material and its general 
concentration from LANDSAT images. To do so one must manipulate the total 
signal from LANDSAT so that only the volume reflectance, py, . is left as 
a residual. Only py relates to the types and concentrations of materials 
in water. Surface effects, bottcm effects, and atmospheric effects are all 
noise sources vdiich must be removed. In order to remove these effects the 
volume reflectance, py, of distilled water or a very clear lake approaching 
distilled water must be determined in the laboratory or by use of the Bendix 
Radiant Power Measuring Instrument (RPMI) (2) . Laboratory values are more 
precise than field values because in the field one must contend with indir- 
ect skylight and wave action which can be eliminated in the laboratory. 

On the LANDSAT image a clear lake which approaches distilled water must be 
used for calibration. With the LANDSAT signal from this lake and the known 
volume reflectance for its clear water it is possible to eliminate the sur- 
face and atmospheric effects and have residual signals which are indicative 
only of the type and concentration of the material in other lakes. 

MADISON AREA LAKES 

Figure 1 shows two LANDSAT images of lakes near Madison, Wisconsin. The 
upper image is taken in early spring shortly after the ice had thawed; 
there was essentially no algae growth in any of the lakes. All lakes 
appear essentially the same brightness on the LANDSAT image. 

The lower image shows the same area in late August when algae and weed growth 
are at a maximum. Those lakes which have an abundance of nutrients (eu- 
trophic lakes) have heavy algae growths. Skylight and sunlight interact 
with this suspended algae and is backscattered to the satellite. The denser 
the algae growth, the more the backscatter and the brighter it appears on 
the LANDSAT image. A clear (oligotrophic lake) such as Devil's Lake does 
not have enough nutrients present to sustain significant algae growth and 
it has approximately the same strength of backscattered signal in August 
as in early spring. The highly eutrophic lakes such as Lake Kegonsa have 
backscattered signals so high in August that they are virtually indisting- 
uishable from the green fields in Figure 1. 

The backscatter from a lake can be sensed in the four bands of the LANDSAT 
satellite. The long -established water quality parameter used to indicate 
backscatter is turbidity. Turbidity is an average value determined across 
the visible spectrum and its units are expressed as FTU's. Figure 2 shows 
the strength of the backscatter signal from Madison area lakes as sensed 
from the LANDSAT satellite on two different days in late summer. These 
signals are plotted against turbidity. It is obvious from Figure 2 that 
for a particular day there is a good correlation between the LANDSAT 
signal and turbidity and eutrophic classification. However, for a differ- 
ent day the height of this correlation curve shifts. This shift is caused 
by different atmospheric parameters from day to day. If one can determine 
the exact location of the curve for a particular day it can be used to map 
turbidity if there are no tannin lakes present nor lakes where bottom ef- 
fects are significant. 

For tannin lakes, the tannic acid in the water creates a dark brown color 
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APRIL 



LATE AUGUST 

FIGURE 1 ERTS (LANDSAT) Images of Lakes near Madison, 

Wise, in early Spring when Algae and Lake Weed 
Growth is a Minimum and in Late August when it 
is Maximum. Band 5 (Red) 
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Satellite Signal, % of Maximum 


Suspended Solids 



.3 ,5 1.0 2.0 3.0 5.0 10 20 40 

TXirbidity, FTU 


• = 28 August 1972 0=5 August 1973 

A = Devil's Lake (clear, oligotrophic lake) 

B = Lake Kegonsa (algae filled and eutrophic) 


Figure 2. Strength of backscatter signal received by Band 5 of ERTS 

(LANDSAT). Madison Area Lakes. Note shift in height of curve 
for different days due to atmospheric change. 
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which absorbs energy. On the curve of Figure 2 such a lake would lie be- 
low the curve for the other lakes. For a lake with strong bottom effects 
the signal from say a sand bottom would create a point which lies far above 
the curve for the other lakes. In order to handle tannin and sand bottom 
lakes the spectral signatures must be analyzed first so the lakes can be 
put into classes of essentially tannin lakes, non-tannin lakes and bottom 
lakes. Within a class the curves, such as those in Figure 2, can be used 
to categorize the lakes as to turbidity and possible eutrophic classifica- 
tion. The lakes in Figure 2 are essentially non-tannin lakes with various 
amounts of algae present and no bottom effects. 

In Figure 2 the shift of the signal -turbidity curve is considerable. There 
is no universal curve for all days. On the other hand, the laboratory curve 
of backscatter versus turbidity is essentially universal for all days. 

Figure 3 shows such a curve which combines laboratory red-light backscatter 
values for 127 samples analyzed over three years at the University of Wisconsin. 
Values for tannin lakes tend to fall slightly below the average curve and 
form a curve of their own, but generally speaking the laboratory backscatter 
versus turbidity curve is a universal curve. The laboratory backscatter 
in Figure 3 is indicated on the y axis by the apparent laboratory reflect- 
ance, AP, where 


AP = Pv PSL 
PPL 

and pp^is the diffuse reflectance of the standard laboratory reflection 
panel. The calculated value of ppL -65 microns was 39%. The average 
approximate value of pcp can be considered about 0.020%. With these values 
one can, with any Pv* obtain the corresponding expected value of turbidity. 
Once turbidity is known it is then possible to map other water quality 
parameters which might correlate to turbidity for that particular type of 
water . 

The correlation of other water quality parameters such as suspended solids 
to turbidity is not universal but varies for different waters. It is possible 
to have a few large particles of dark material which have a certain weight 
which scatter back considerably less energy than say a large number of 
exceedingly fine white particles of the same total weight. For the same 
sized and shaped particles of a certain material there is a correlation 
between weight of suspended material and turbidity, but one must point out 
that this correlation will not necessarily hold for another material in a 
different type water. 

When the material causing the turbidity or backscatter is say algae then 
there may also be a correlation between backscatter (satellite signal 
strength) and chlorophyll. But try the same correlation on an inorganic 
red clay and the backscatter-chlorophyll correlation breaks down. If the 
suspended material causing turbidity is municipal or industrial sewage 
which upon decay uses oxygen then for this situation there may be a correl- 
ation between backscatter (remote sensing signal strength) and Biological 
Oxygen Demand (BOD) , or the dissolved oxygen which is not yet used up (3) . 
However, try the same correlation on inorganic red clay in water and the 
correlation breaks down. For all v;aters however, the correlation of back- 
scattcr (volume reflectance, ) and turbidity holds as showT. ir. Figure 5. 

For Red Energy (.65 microns) this correlation is; T = 5.21(AP)2.00 
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Figure 3. Laboratory backscatter expressed as apparent laboratory reflectance, AP, plotted 
against turbidity for 127 different lake samples collected over three years. Th 
curve is for red light (0.65 microns). The best fit equation for this data is 



For this investigation, the volume reflectance values, p , for a water are 
obtained for different wavelengths, some of which corrcs'|Sond to the LANDSAT 
bands. The resulting spectral signature allows one to jiut lakes into general 
categories such as tannin or non-tannin waters. Then within a category 
a curve such as in Figure 3 can be used to categorize the lakes as to re- 
lative turbidity and the amount of suspended material such as algae v,hich 
causes this turbidity. 


LABORATORY ANALYSIS 

The volume ref lectance,p.^. , can be obtained by the scheme shown in Figure 4. 
In a-darkened laboratory, a laboratory lamp provides Tight onto a standard 
BaSO^ diffuse reflection panel and onto the water sample. The water sample 
is in a tube .62 meters long which has a black bottom. Shorter tubes will 
cause overriding bottom signals which make the results confusing or meaning- 
less as will a bottom which is too reflective (4) (5). The sides of the 
tube are lined with diffusely reflecting white chronoflex which returns 
side-scattered energy simulating adjacent water volumes in the field. 

The radiance of the lamp is in watts/cm2 s. As viewed from the level 
of the panel or water sample the solid angle of the lamp is steradians. 
The total irradiance that reaches the level of the panel at right angles 
to the rays is . The total irradiance available on the panel per unit 
area is 


= Lj^iJ/j^cose (Vfltts/cm2) 

where 0 is the angle between the lamp and the vertical as shown in Figure 
4. The indirect illumination from the ceiling L^. (watts/cm^ s) is zero 
in the perfect laboratory setup. 

Lambert’s Law states that P, the diffuse radiance returning into space from 
the panel, is 


P = equation 1 

•n 

where pp^ is the diffuse reflectance of the panel. 

Since the dust, foam and other impurities on the water surface also behave 
as diffuse reflectors their signal, S^. also follows Lambert's Law 

ir 

where P$i is the diffuse reflectance of the impurities on the water surface 
in the laboratory. The suspended particles in the water volume are also 
diffuse reflectors and cause the radiance from the water volume V. 

IT 

where, P.^. = the diffuse reflectance of the particles in the water volume. 
This Py is the all-in^jortant factor which is indicative of water quality. 
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Figure 4. Direct Energy Relationships in the Laboratory 



The signal G does not follow the laws of diffuse reflection. When e = r, 
as in Figure 4, the signal G is specular reflection of the lamp's radiance 
from the air-water interface. 


G = (Natts/cirrs) 

tdiere 4> = the surface specular reflection of a water-air interface. 


Ha 


2 

^ _ (^w where 

("w "a) 2 

index of refraction of water = 1.333, and 
index of refraction of air = 1.000. 


0 


(1.333 - 1.000)^ _ 
(1 .333 + l.CC0)2 " 


This value holds fairly well out to angles of 6 and r of about 40 degrees 
of the vertical. Most remote sensing operations fall within this range. 


The total signal from the water as seen by the sensor is w 


W = 
W = 


V + s. 




» = Cp^ - Pjl) I!k 

IT 


equation 2 


The ratio of this signal W to the signal from the panel is called apparent 
laboratory reflectance AP. 

equation 3 


The PpL reflectance for BaSO^, for 0.65 microns by use of equation 1, was 
calculated to be 39% ^ 1.5%. Therefore AP = (p + c )/.39. Figure 3 is 
a plot of this AP for waters of different turbidities. 

Let subscript "i" indicate values for distilled water and subscript ." 2 ” 
indicate values from water sample #2 (a more turbid water) . Ke can define 
A ?2 ~ APj as the ''Laboratory Difference, D 2 " between the laboratory 
apparent reflectances of samples #2 and #1 (distilled water). 


AP = — = 

P ' 




(p ^ 


D2 = AP2 - APj 

D 2 = ^ ^SL 

PpL 




equation 4 


_ — .^t^rirroTT .TTY 


ri"l 


10 



From Figure 3 vie see that for distilled water the average value of APj^ 

/J>, = ^vl ^SL = 0.18% 

PpL 

A usable value for Pc^ is .020%. This leaves a value for Pyj = .18% (.39) -.020% 
P yj - .050%. 

From equation 4 if we have 02; we can calculate Py2 by 

%2 ^ ^2^ri ^vl 
Pv2 =02-39 + .05% 

The value 

D 2 = ^v2 ~ Pyl 
^PL 

is the important factor which will be further used to correlate between 
laboratory values and satellite values. 

BOAT ANALYSIS 

When the water being analyzed is out of doors as in Figure 5 we still have 
the direct radiance of the sun Lg, which behaves analagously to the labora- 
tory lamp in Figure 4. l\'e also have L^, the average indirect radiance from 
the skylight. According to Lambert's Law the total irradiance effecting a 
horizontal surface due to this skylight is 

He » 

The sun's direct contribution is = Ls^cosO. Let the total irradiance 
be h; 

I t . I 

+ 11 

o S c 

t I f 

= L,i|cos€’ + L It 
os c 

The diffuse signal from a field reflection panel, P, would be 

P = 


idio'e Pp = the reflectance of this field panel. 

Also following Lambert's Law for diffuse reflection, S^, the diffuse reflection 
caused by impurities on the water surface is: 

S* = 

s 

It 


where Pj = the field diffuse reflection for foam, leaves, and other dirt 
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Figure 5. Direct and Indirect Energy Relationships in the Field. 


on the water surface. Also V*, the radiance from the water volume is 

y ’ _ 

■n 

The specular reflection of the sun from the water surface (sun glitter) is 

G' = Ls« = .020 Lg 

In a similar ^manner the specular reflection of the skylight from the water 
surface is S^. 

Sc = Let = .020 Lc 

The total signal from the water is W* . 

W' = V' + Sg + Sc 

' D F' o H' • 

K = o + o + .020 L 

c 

ir Tt 

k’ 

K' ^ — + .020 l' 

n c 

The signal from the field panel is 



If one were determining the volume reflectance of a particular lake in the 
field with the Bendix RPMI ground truth instrument this can most easily 
be done by pointing the RPMI sensor telescope down through the surface of 
the water as in position "d" in Figure 5. Assuming no significant bottom 
signal the sensor at "d" only reads the radiance from the water volume v'. 


Hq can be obtained directly by use of the RPMI or a reading can be taken 
on the panel. 



— — and 

IT Pp ■ 

p = V Pr 

, V 

Thus the Py can also be obtained in the field for any water san 5 )le with the 
RPMI. However, due to variations in skylight and wave action in the field 
the laboratory determination of is a magnitude more precise than field 
determinations. 
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AIRBORNE ANALYSIS 


When the signal from the water level such as Sg in Figure 5 passes upward 
through the atmosphere it is attenuated by the atmospheric transmittance x ; 
^ = SgX. Also the atmospheric backscatter, LA, is added to the signal. 
Therefore, the total signal from the water as sensed at "c" in Figure 5, is 


W" = 


w 


w" = 


(v' + Sg + S^)t + LA 

pH* + pH' + .020 C 

VO so ^ 


* LA 


^ ^ H 

(Py Ps^ 

IT 


+ .020 L T + LA 

c 


If we let subscripts i and 2 represent distilled or very clear water and 
a turbid lake (#2) respectively, then 

I 

Wj =(Pyj + Pg)!!^-^ *020 L^T + u 
ir 

W '2 =(Py2 Ps^— -020 L^t + lA 

TT 

Also if we define the satellite residual R '2 to be the difference in signals 
which is caused only by the material in water #2, then 


R 2 = W 2 - Wj = (p, 


v2 




Again assuming the volume reflectance of distilled water, Pyj, to be known; 


v2 


= r: 


Ho- 


vl 


equation 5 


With this Py2 one 
What is needed in 
Then this reading 
Wi. 


can enter the curve on Figure 3 and obtain turbidity . 
the above situation is w'l, a reading on a very clear lake, 
can be subtracted from the reading from any other lake 


K'l = Wi - wi’ 


and 


P 


vi 




+ 


%1 


If is not known the answers are in relative terms only, 
of a turbid lake (sample #2) is known then from equation 5, 


If the Py 2 


IT 

This H^term can then be used 
of any lake. 


IT = Py2 



to solve for the absolute volume reflectance 
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The results in this paper are for relative values of Pyi only. The exact 
mathematical relationships were not fully understood at the time of the 
overflights in order to obtain the necessary simultaneous ground truth 
for sample #2. 


SATELLITE ANALYSIS 

On the LANDSAT image assume we can locate a very clear lake; ”#1”, which 
approaches distilled water in purity. This lake must be deep enough so 
that no bottom is showing. Aerial observation should be made to assure that 
there are no bottom effects showing in this test lake. 

Let Wj equal the satellite raw reading on the clear lake. Let W^' equal 
the raw reading on Lake #2, the lake in question. The raw readings for 
Band 4 and Band S are very high due to atmospheric effects and the different 
ial water signals are perhaps 1% of these raw readings. The raw readings 
are not very meaningful . lbv;ever, if we subtract Vi’{ from U '2 the residual 
is due only to the material in Lake #2. 


Rj = k” - w;' = - ‘>vP«o i- 


The four bands of LANDSAT cover as follows: 
Band 4: 0.50 to 0.60 microns 

Band 5: 0.60 to 0.70 microns 

Band 6: 0.70 to 0.80 microns 


Band 7: 0.80 to 1.10 microns 

Let us plot the total signal from Bands 4, 5, 6, and 7, as centering on 
wavelengths 0.55, 0.65, 0.75, and 0.95 respectively. Figure 6 shows such 
a plot of the satellite residual, R^, for these four LANDSAT bands. Three 
types of waters; clear, tannin, and algal, are readily differentiated. 
Therefore, it is possible to separate lakes into these classes by LANDSAT 
signals alone. Also, within a particular class such as algal, one can 
categorize the lake as to how much algae is present. Figure 7 shows 
satellite residual signals of lakes which grade from clear to heavy algal. 

The shape of the satellite residual curv'es of Figure 7 duplicate well the 
shape of the laboratory reflectance difference curves in Figure 8 for 
the same types of water. In Figure 7, the satellite residual 


R. = (p . 

1 ’“VI 


P -I )H - 

Vl-' O TT 


In Figure 8, the laboratory reflectance differences are: 




P 1 )-i- = (P • ■ P 1 )-L 
vl''~ '■VI vl'' — 


PpL •" .39 

Therefore, the values on the two figures vary only by tlie factor: 


p** ^ T 1 ® T 

= H - T -i- = .39 H 1 


cr 


o r 


.39 


O TT 
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Satellite Residual Signal, R^, LANDSAT Sensor Counts 



.55 .65 .75 .95 

Wave Length, microns 


Figure 6. Satellite Residual Signal, R^, for three v;ater types. 

R- . k’ - w" . Cp„. - 


16 


Satellite Residual Signal *'Ri" LANDSAT Sensor Count 



T = turbidity 

Figure 7. Satellite Residual Signal "Ri" for clear-type water 
lakes with various amounts of algae present. 


R. = (p . - p -)H f 
1 '■VI ^vl o — 

IT 
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Laboratory Reflectance Difference, D., Percent 


e 
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0.1 
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- 0.1 



Figure 






The gradation between heavy algal and clear water lakes are equally apparent 
in both Figures 7 and 8. Exact corresponding turbidities arc also shown 
in Figure 8. Three approximate turbidities are shown in Figure 7. On 
Figure 8 one can also see that it is possible to tell the differences be- 
tween types of algae. Blue-green and green algae types are indicated. 

Green algae are usually not the nuisance types but blue-green algaes are^ 
and are often associated with eutrophic lakes which can become low in 
dissolved oxygen. 

Figure 9 shows that although the tannin type lakes have a different shaped 
characteristic curve, the more turbid tannin lakes have a higher signal and 
follow the same general pattern as \dth the algal lakes. Figure 9 also 
shows the relative ranges of the satellite and laboratory data. 

OPTIMUM TIME OF YEAR FOR EUTROPHIC CLASSIFICATION 

Figure 10 shows that in an individual lake the volume reflectance can change 
markedly from a characteristic clear water lake in early spring to a very 
algal lake in late summer and back to a clear water type lake in late fall. 

This is also indicated by Figure 1. The maximum algal and weed growth is 
in late August. This is also the time of minimum dissolved oxygen (see 
Figure 10) and also the time of maximum water temperature and when the nitrates 
become completely tied up in some lakes. From many indicators, late August 
is the optimum time for eutrophic classification of lakes. 

BOTTOM PROBLEMS 

The bottom effect problem is a difficult noise factor in satellite classi- 
fication of lakes. Bottom signals can show in any of the lake types. The 
type of bottom (dark mud, light sand, or green weeds), also will give a 
different characteristic modification to signals from each of these water 
types. Also the depth to bottom effects the strength of the signal and 
also its spectral distribution. 

It is anticipated that the bottom effect problem can be completely isolated 
by analyzing LANDSAT data from spring imagery when the lakes are clear and 
free of algae and the bottoms are very apparent. Where bottom signals are 
strong they are a characteristic signal and can be classified as another lake 
type. Figure 11 shows sand bottom and weed bottom effects on the satellite 
signals from a tannin type lake. 

COLOR CATEGORIZED IMAGES PRODUCED BY THE 
BENDIX CORPORATION M-DAS EQUIPMENT 

By training the Bendix M-DAS computer for each type of lake of interest 
the machine recognized all such types in the scene and displayed desired 
types as a color categorized image. It must be emphasized that good ground 
truth in the form of some water samples and aerial observations of the test 
lakes are essential for the training of the M-DAS equipment. The aerial 
observations are especially essential in locating the troublesome bottom 
effects which might show up on a training lake chosen without aerial observation. 
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Laboratory Difference, D^, percent 



Laboratory Reflectance, AP, 





Wavelength, microns 



DO = Dissolved oxygen 


Figure 10. Time of year effect on Laboratory Reflectance of Madison 
Area Lake. Maximum Reflectance occurs at the same time 
as minimum dissolved oxygen, which is late August. 
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Figure 12 Initial Color Categorized Image Produced by Bendix Corp. 

Showing Lake Categorization. (LANDSAT Scene 1756-16061 
of 18 August 1974) 
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Figure 12 shows the test color categorized image for lake categorization 
produced by the M-DAS equipment. The results look good and are being field 
checked in August, 1975. Further anticipated work includes combining spring 
imagery vdth summer imagery to completely solve the bottom effect problem. 
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